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ABSTRACT 

Aims. The postburst object of a GRB is likely to be a highly magnetized, rapidly rotating compact object (e.g., a millisecond 
magnetar), which could produce an ultrarelativistic electron-positron-pair wind. The interaction of such a wind with an 
outwardly expanding fireball ejected during the burst leads to a relativistic wind bubble (RWB). We investigate the properties 
of RWBs and use this model to explain the shallow decay phase of the early X-ray afterglows observed by Swift. 
l\/lethods. We numerically calculate the dynamics and radiative properties of RWBs. 

Results. We find that RWBs can fall into two types: forward-shock-dominated and reverse-shock-dominated bubbles. Their 
radiation during a period of ~ lO'^ — 10^ seconds is dominated by the shocked medium and the shocked wind, respectively, 
based on different magnetic energy fractions of the shocked materials. For both types, the resulting light curves always have 
a shallow decay phase, as discovered by Swift. In addition, we provide an example fit to the X-ray afterglows of GRB 060813 
and GRB 060814 and show that they could be produced by forward-shock-dominated and reverse-shock-dominated bubbles, 
respectively. This implies that, for some early afterglows (e.g., GRB 060814), the long-lasting reverse shock emission is strong 
enough to explain their shallow decay phase. 
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1. Introduction 

One of the most puzzling features of the early X-ray af- 
terglow light curves of the gamma-ray bursts (GRBs) dis- 
covered by Swift (Gehrels et al. 2004) is the existence of a 
flattening segment (temporal indices a ~ [0, —0.8]), which 
lasts from a few hundred seconds to a few hours (Campana 
et al. 2005; Vaughan et al. 2005; Cusumano et al. 2005; 
Nousek et al. 2006; O'Brien et al. 2006; de Pasquale et al. 
2006; Willingale et al. 2006). This feature has been widely 
understood as due to a long-lasting energy injection. Two 
kinds of energy injection have been proposed. One kind is 
the so-called refreshed shock scenario with a smooth dis- 
tribution of the Lorentz factors of the ejected shells (Rees 
& Meszaros 1998). The other kind, in focus here, involves 
a central engine activity extending over a long period (Dai 
& Lu 1998a, 1998b; Zhang & Meszaros 2001; Wang & Dai 
2001; Dai 2004, D04 hereafter; Zhang et al. 2006; Fan & 
Xu 2006). 

The popular models for the origin of long and short 
GRBs are the collapse of a massive star and merger of 
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a compact binary, respectively (for recent reviews see 
Woosley & Bloom 2006; Nakar 2007). These models pre- 
dict that, after a GRB, the remaining compact object 
seems to be a millisecond-period pulsar or a rapidly- 
rotating black hole. If the pulsar is strongly magnetized 
(e.g., a magnetar) or the black hole has an accretion disk 
lasting for a long time, these compact objects will con- 
tinuously release their rotational energy through some 
magnetically-driven processes to produce an energy out- 
flow. As this magnetically-driven outflow catches up to 
and then interacts with the relativistic fireball ejected dur- 
ing the GRB, the fireball's energy will increase. 

Based on this consideration, Dai & Lu (1998a, 1998b) 
and Zhang & Meszaros (2001) propose an energy-injection 
model for GRB afterglows with an assumption that the 
energy outflow is purely composed of low-frequency elec- 
tromagnetic (EM) waves radiated by a postburst magne- 
tar. Using this pure EM energy-injection (PEMI hereafter) 
model, Fan & Xu (2006) successfully explain the shal- 
low decay phase of the X-ray afterglow of GRB 051221a. 
However, the numerical calculations (Wang & Dai 2001; 
Fan & Xu 2006) also show that an actual flattening seg- 
ment of a light curve is steeper than an analytical estima- 
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tion when the dynamics is numerically described and th( 
equal-arrival surface effect is considered. Thus, it couk 
be difficult to use the PEMI model to explain some GRE 
X-ray afterglows with a very flat plateau, such as GRE 
060814. 

On the other hand, the PEMI model does not con- 
sider possible evolution of the energy outflow with ra 
dius. Because the fluctuating component of the magnetic 
field in the outflow can in principle be dissipated by mag 
netic reconnection and used to accelerate an associatec 
electron-positron plasma, the outflow should eventually 
become a kinetic-energy flow carried by the acceleratec 
e^ pairs, even though it is dominated by the EM en- 
ergy at small radii (Coroniti 1990; Michel 1994; Kirk & 
Skjaeraasen 2003). This transformation of the outflow is 
estimated to occur around the radius ~ lO^'rL — lO^rL, thai 
is, below the deceleration radius of the fireball ejected dur- 
ing the burst, where tl (~ lO^cm for millisecond magne- 
tars) is the light-cylinder radius of the pulsar. With these 
arguments, D04 suggests that it is likely to be an ultra- 
relativistic e^-pair wind (but not pure EM waves) that 
interacts with the fireball to influence the GRB afterglow. 
As a result of the interaction, a rclativistic wind bubble 
(RWB hereafter) is produced, which is a relativistic ver- 
sion of the Crab nebula. 

In this paper, we calculate numerically the dynamic 
evolution and the corresponding radiation of a RWB 
driven by a millisecond magnetar by considering the 
inversc-Compton scattering effect and the equal-arrival 
surface effect of the RWB. Because of these effects, our 
numerical results are different from the analytical ones of 
D04. In particular, we find that RWBs fall into two types, 
which are dependent of the magnetic energy fractions of 
shocked materials. As an example, we fit the observed X- 
ray afterglows of GRB 060813 and GRB 060814, which 
are found to belong to two types of RWBs. We emphasize 
that, for some early afterglows (e.g., GRB 060814), the 
long-lasting reverse shock emission is strong enough to 
explain their shallow decay phase, which cannot be simu- 
lated by using the PEMI model. 

2. Dynamics of a RWB 

Some energy-source models (for a review see Zhang & 
Meszfiros 2004) suggest that the central engine of a GRB 
is a millisecond magnetar. In such models, a GRB itself 
may be duo to neutrino and/or magnetic processes of a 
rapidly-rotating, strongly-magnetized pulsar, which may 
eject a few shells. Collisions between the shells may pro- 
duce internal shocks, which give rise to a few pulses of 
an observed GRB during the prompt emission. After the 
GRB, this pulsar will be braked by magnetic dipole radia- 
tion. As a result, the release of the stellar rotational energy 
(at a rate of i^rot) drives an outflow, whose luminosity is 
estimated by (Shapiro & Teukolsky 1983) 

Lw = —Erot 

^ 4.0 X 10^'Bl,,Rl,P,-^^ (l + ^)"' erg s'S 




Fig. 1. A schematic cartoon of an RWB. The meaning of 
the four regions (1-4) is explained in the text. In the in- 
ner (a) part of region 4, the energy outflow is dominated 
by low-frequency EM waves, in the outer (b) part of re- 
gion 4 by a kinetic-energy flow carried by ultrarelativistic 
electron-positron pairs. 



where Sj__i4 = Bs^i4sin^, 9 is the angle between the 
magnetic and rotational axes, and -Bs,i4 = -Bs/IO^^G, 
-Rs,6 = i?s/10^cm, and Po,ms = Po/lms are the surface 
magnetic field strength, the radius, and the initial pe- 
riod of the magnetar, respectively. The time t is mea- 
sured in the observer's frame. The characteristic spin- 
down time is T^d ~ 5.0 x 10^^(1 + z)Bl^^Ji5R-lPl^^ s, 
where 145 = I/IO^^ g cm^ is the stellar moment of inertia 
and z the cosmological redshift. In the following calcu- 
lations, the typical values of i?s,6> -Pojins) and 745 are all 
taken to be unity. 

As the outflow propagates outward, the energy dissipa- 
tion of the EM-wave component significantly accelerates 
plenty of associated pairs. Eventually, the outflow be- 
comes an ultrarelativistic e^-pair wind around some cer- 
tain radius that is below the deceleration radius of the 
shells ejected during the burst. The wind's bulk Lorentz 
factor is 7w ^ 10^ — 10^, as argued by Atoyan (1999) for 
the Crab pulsar. In the following calculations 7w = 10''. 
However, we find that a higher value of 7w does not change 
our results significantly. As a result of interaction be- 
tween the wind and the medium, an RWB should include 
two shocks: a reverse shock that propagates into the cold 
wind and a forward shock that propagates into the am- 
bient medium. Here, for simplicity, we assume that two 
initially-forming forward shocks during the interaction of 
the GRB ejecta, both with the medium and with the wind, 
have eventually merged to one forward shock. Therefore, 
as illustrated in Fig. 1, there are four regions separated 
in the bubble by the shocks: (1) the unshocked medium, 
(2) the forward-shocked medium, (3) the reverse-shocked 
wind gas, and (4) the unshocked cold wind, where regions 
2 and 3 are separated by a contact discontinuity. 

We denote some quantities of region i as follows: is 
the particle (proton or electron) number density and Pi 
the pressure measured in its own rest frame, and 7^ and 
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/3j = (1 — 7j~^)^^^ are the bulk Lorentz factor and cor- 
responding velocity measured in the local medium's rest 
frame, respectively. The total kinetic energy of region 2 
is Ek,2 = 72(mej + msw)c^ + 72(72 - l)mswC^, where rriej 
is the rest mass of the initial GRB ejecta and nisw is the 
rest mass of the swept-up medium. Energy conservation 
requires that any increase in Ek,2 should be equal to work 
done by region 3: dEK.2 = 5W = AirB^P^dR, where R is 
the radius of the bubble in the thin shell approximation. 
Then, we can obtain 



dR ~ 



(7l 



l)nimp 



mej + 272msw 



(2) 



On the other hand, the dynamic evolution of region 3 can 

be determined by the relationship between the Lorentz 
factors of the two sides of the contact discontinuity surface 
according to Blandford & McKee (1976) 



73 = 72X 



-1/2 



where the similarity variable x is (D04) 



X = 



1 67rn 1 mp c-^ 7I i?2 



-12/29 



(3) 



(4) 



where rUp is the proton rest mass. Using Eqs. 2, 3, 4, P3 = 
Lw/127ri?2-y2(, (D04), and the relations oi dt = (l + z)(l- 
f33)dR/l33C for region 3 and dt = {1 + z){l - f32)dR/l32C 
for region 2, we can get the dynamic evolution of the 
RWB before the characteristic time Tgd- Simultaneously, 
the masses of the shocked medium and shocked wind gas 
are calculated respectively by 



drus 



dR 
dm^ 



= 47ri?^nimp, 



= 47ri?^(/^4 - /3Rs)74"4TOe, 



(5) 



(6) 



where rric is the electron rest mass and fSj^g = (73^3/33 — 
74714/34)/ (73713 — 74/14) is the velocity of the reverse shock 
measured in the local medium's rest frame. However, when 
t > Tgd, the reverse shock is regarded as terminativc 
(D04) and thus Eq. 3 should be replaced by 73 ex R^''/'^ 
(Kobayashi & Sari 2000; Kobayashi 2000), as well as 
ns oc i?-i3/2^ P3 cx i?^26/3^ ^ 0. 

Figure 2 shows the evolution of the bulk Lorentz fac- 
tors of regions 2 and 3 by taking the initial isotropic ki- 
netic energy Ek.2.o = lO^'^ ergs and the initial Lorentz 
factor 72,0 = 150 of region 2 at the deceleration ra- 
dius, and ni = 1 cm~^. The shape of the curves is ob- 
viously consistent with the analysis in D04, who points 
out that the evolution can be divided into three stages: 

(I) 72 oc i-3/8,73 OC t-39/136. (n) 72 ~ 73 OC i'V^; 

(III) 72 CX i~^/*,73 cx t^'^/^^'. Compared with the case 
without wind injection, a significant change in the evo- 
lution of 72 is the existence of stage II, which is quite 
similar to the result of the PEMI model (Wang & Dai 
2001). This similarity can be easily understood since the 
term ^TrR^P^dR = |(1 -|- z)~^L^dt, which represents the 
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Fig. 2. The evolution of 72 and 73 in the RWB model 
with B_L44 = 2 and z = 1. The dotted curve represents 
the bulk Lorentz factor of the shocked medium without 
wind injection. The initial value of 72 is taken to be 150. 
This order of magnitude has been suggested by some ob- 
servations (e.g., Molinari et al. 2006) 



work done by region 3 to region 2 in Eq. 2, is approx- 
imately equal to (strictly, a factor 2/3 times) the term 
z)~^ L^dt in the dynamic equation of the PEMI model 
for the same magnetar (for details, see Dai & Lu 1998a, 
1998b; Wang & Dai 2001; Fan & Xu 2006). Therefore, 
we cannot expect to distinguish the RWB model from the 
PEMI model only by observing the radiation from the 
forward-shocked medium. A possible difference between 
these two energy-injection models is induced by reverse- 
shocked pairs, the radiation of which arises from the re- 
maining one- third energy release of the magnetar. 

3. Radiation and example fit to X-ray afterglows 

Both the forward shock and the reverse shock heat cold 

materials to a higher temperature, generate random mag- 
netic fields, and accelerate protons and electrons. Since 
the microphysical processes have been unclear so far, the 
electron energy density and the magnetic energy density 
are parameterized as usual. We assume that for region 3 
the electron and magnetic field energy densities are frac- 
tions 60. R and eb.r of the total energy density behind the 
reverse shock (ee,R -|- eB,R = 1), and for region 2, frac- 
tions Ce.F and CB.F of the total energy density behind the 
forward shock, where ec,F + (b.f < 1 and ec,F ^ ^ee.F 
according to Medvedev (2006). It is natural to think that 
eb.r 7^ cb.f and ee,R ^ ee,F, as suggested in some studies 
(Fan et al 2002; Coburn & Boggs 2003; Zhang, Kobayashi 
& Meszaros 2003; Kumar & Panaitescu 2003). We also 
assume that the spectral indices of the electron energy 
distribution are p2 and for regions 2 and 3, respec- 
tively. Furthermore, we assume here that P2 ~ P3 = P- By 
giving a set of three free parameters {p, cb.f and cb.r). 
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Fig. 3. Parameter space of {p, cb.f, cb.r) with B±s4 = 
2 and z = 1. The hnes represent the critical values 
(eB,R/eB,F)c with different p, below which ( in the shaded 
area for p = 2.3) the corresponding reverse shock emis- 
sion is insignificant for the total fltix of the RWB. The 
star symbols label the parameter sets of two GRBs. 



we can fix the cooling Lorentz factors 7e,c arising from 
both synchrotron radiation and inverse-Compton scatter- 
ing, and the niinimuni Lorentz fcictors '7c.min 

in the elec- 
tron energy distributions of regions 2 and 3. We also con- 
sider the "equal-arrival surface" effect of the RWB emis- 
sion. The related formulas for these calculations were pre- 
sented in Sect. 3 of Huang et al. (2000). We ignore the 
self- absorption effect for the X-ray emission we are inter- 
ested in. 

Before light curves are exhibited, it should be pointed 
out that the ratio of the flux contributed by regions 2 or 
3 to the total flux is sensitive to the ratio of eB,R/eB,F- 
As shown in Fig. 3, if a dot representing a parameter set 
{p, eB,F, cBjr) is in the shaded area (for p = 2.3), the 
peak of the light curve of the reverse shock emission is 
below the light curve of the forward shock emission (e.g., 
see the upper panel of Fig. 4). The higher the value of p, 
the easier the occurrence of this case. Because the forward 
shock dominates the radiation of the RWB, as discussed 
above, the combined light curve should be similar to the 
result of the PEMI model for an identical magnetar. In 
other words, all the GRB X-ray afterglows that can be 
fitted in the PEMI model must be explained by the RWB 
model with an appropriate parameter set. Conversely, if 
eB,R/eB,F is large enough (> (eB,R/eB,F)c)) tiie radiation 
from region 3 becomes quite important, especially dur- 
ing a period of ~ 10^ — 10^ s. Therefore, we can divide 
the RWBs into two types: forward-shock-dominated and 
reverse-shock-dominated bubbles. 

However, because of the existence of stage II in the dy- 
namic evolution, the afterglows produced by two types of 
RWBs always have a shallow decay phase. This feature of 



RWBs is very consistent with a lot of afterglows observed 
by Swift. For example, two X-ray afterglows of GRB 
060813 (upper panel) and GRB 060814 (lower panel) are 
fitted by the RWB model as shown in Fig. 4, and the corre- 
sponding parameter sets {p, ee.F, ^b.r) are shown in Fig. 

3. It can clearly be suggested that the afterglows of GRB 

060813 and GRB 060814 could be produced by forward- 
shock-dominated and reverse-shock-dominated RWBs, re- 
spectively. As discussed above, the prompt emission of a 
GRB may result from internal shocks. If the high-latitude 
emission of the last internal shocks is weaker than the early 
afterglow emission, one can observe the shallow decay of 
an early afterglow, as in GRB 060813. Conversely, an early 
steep decay should be observed (Kumar & Panaitescu 
2000), as in GRB 060814. For the latter burst, we should 
thus introduce this curvature effect to explain the earlier 
steep decay, which is fitted roughly by the following equa- 
tion (as shown in the inserting panel of Fig. 4, also see 
Zhang et al. 2006; Liang et al. 2006): 

F^{t) = A(^\ +Fx,rwb, (7) 

where A is the normalization parameters, to the time zero 
point of the last prompt emission pulse, and /3 the spectral 
index that is equal to 0.67 (Moretti, Guidorzi & Romano 
2006). Please note that we take the time zero point of 
the afterglow as the GRB triggering time as proved by 
Kobayashi & Zhang (2006) and Liang et al. (2006). 

4. Summary 

In this paper, we have numerically calculated the dynamic 
evolution and radiation of RWBs. Two motivations impel 
us to consider this RWB model rather than the PEMI 
model: first, the PEMI model cannot explain the very 
flat plateau in light curves of some X-ray afterglows; sec- 
ond, a pulsar-driven energy outflow, which is dominated 
by Poynting flux at smaller radii, could evolve into an ul- 
trarelativistic electron-positron-pair wind at larger radii. 
The most significant feature of the light curves of the X- 
ray afterglows from the RWBs is a flattening segment oc- 
curring during a period of ~ 10^ — 10^ s, which is consis- 
tent with the observed shallow decay phase. Our exam- 
ple flts to the X-ray afterglows of GRB 060813 and GRB 

060814 indicate that they could be produced by forward- 
shock-dominated and reverse-shock-dominated RWBs, re- 
spectively. This suggests that the central engines of these 
two GRBs could be millisecond magnetars. Moreover, the 
example fits also show that, besides the X-ray afterglows 
that can be fitted in the PEMI model (e.g., GRB 060813), 
the RWB model can also explain some other afterglows 
with a very flat plateau (e.g., GRB 060814) that create 
difficulties for the PEMI model. 

In addition, we would like to point out that the 
magnetar-driven RWB model is one possibility. As pointed 
out in the introduction and in D04, RWBs could also 
be produced, in principle, by central black holes that 
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Fig. 4. Fits of GRB 060813 with B^^u = 6.5, p = 
2.5,eB,F = 0.01, ee.R = 0.0001, z = 0.6, and GRB 060814 
with 5^,14 = 2.5,p= 2.2,eB,F = 0.002, cb.r = 0.0004, z = 
0.6 using the RWB model. The half-opening angles of 
jets are 6j = 0.1 and 9j = 0.15, respectively. The solid 
lines correspond to total fluxes, while the dotted and 
dashed lines represent the fluxes determined by the for- 
ward and reverse shocks, respectively. The early steep de- 
cay segment of GRB 060814 is fitted with A = 9.85 x 
10~^ erg s~^cm~^ and to — 75 s in the insert. 



accrete circumburst materials after the GRB. Moreover, 
some other models have also been proposed to explain 
the observed shallow decay phase, e.g., off-beam jets, 
two-component jets, and varying microphysics parameters 
(Eichler & Granot 2006; Granot et al. 2006; Panaitescu et 
al. 2006; also see Zhang 2007 for a recent review). It is a 
demanding task to distinguish between the above models 
by using multiwavelength observations. 
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